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In order to obtain information about the internal structure of the preformed pairs in the pseudogap 
state of high Tc superconductors, we calculate the propagator of a singlet pair with center of mass 
coordinate r, and relative distance p, by solving the Bethe-Salpeter equation, representing the sum 
over repeated two-particle scattering events due to a distance dependent attraction. We define then 
a "pair structure function" gp{P,p) that depends on the internal distance p between the partners 
and on the momentum P of the pair. We calculate this function both for a local potential and 
s— wave symmetry of the order parameter and for a separable potential and d— wave symmetry of 
the order parameter. The influence of the center of mass momentum, strenght of the interaction, 
temperature, density of particles and of the pseudogap in the one-electron spectrum is studied for 
both cases. 



PACS numbers: 74.20.Fg, 74.20.Rp 
I. INTRODUCTION 

High-temperature cuprate superconductors show re- 
markable deviations from Fermi liquid behaviour in their 
normal state, in particular in the underdoped region. 
One of the most striking features is the opening of a 
pseudogap above the superconducting critical temper- 
ature Tc and below another temperature, T*, that in- 
creases when the doping is reduced^. This pseudogap 
appears to have the same angular dependence and mag- 
nitude as the superconducting gap below Tc, with which 
it seems to merge at the phase transition. In the pseu- 
dogap region, thermodynamic quantities and transport 
coefficients also deviate from Fermi liquid behaviour^. 
To explain these anomalies, different scenarios can be 
invoked. One of them is based on the formation in the 
pseudogap region of incoherent Cooper pairs'^, also called 
"preformed" pairs^. At a lower temperature, phase co- 
herence is established among these pairs leading to the 
superconducting transition. In this respect, the super- 
conducting transition can be regarded almost as a Bose 
condensation of preformed Cooper pairs^, while the BCS 
behaviour, in which the formation of the Cooper pairs 
and their condensation occur at the same temperature, 
is recovered as one approaches the overdoped region. The 
fact that the low temperature coherence length is much 
shorter in the underdoped regime is consistent with the 
picture of preformed pairs well localized in space. ^ In this 
framework, it is interesting to investigate the structure 
of these preformed pairs. 

In the superconducting state, the Cooper pairs can be 
characterized by the so-called "pair-size" £,pair^ , which 
is obtained from the pair correlation function for parti- 
cles paired with opposite spins. In the zero temperature 
limit, this reduces to the modulus square of the anoma- 
lous Green function. This definition can be generalized 
to be valid in the pseudogap region, where the anomalous 



Green function is zero. 

We base our calculations on the two-dimensionale at- 
tractive Hubbard model with an attraction acting either 
on electrons on the same lattice site or on neighbouring 
sites. Analytic and numerical work has shown that such 
a model can decribe the formation of electronic pairs of 
s— wave and d— wave, respectively, and the transition 
to the superconducting state. Going beyond the weak 
coupling limit, the formation af a pseudogap above the 
transition temperature is found. In order to describe the 
internal structure of preformed pairs in underdoped high- 
temperature superconductors we study the Green func- 
tion for the operator which creates a pair of particles 
whose partners are at a certain relative distance, and 
study how this internal distance varies as a function of 
temperature, interaction strength and density of parti- 
cles. To this end, we will consider a Bethe-Salpeter ap- 
proach and introduce a T— matrix formalism, which takes 
into account the relevant interaction channel of repeated 
two-particle scattering. 

In the next section we give an overview of our theoret- 
ical approach. In the case of a dilute system of charge 
carriers, which is of interest to us, the simplest approxi- 
mation to the Bethe-Salpeter equation allows to separate 
the center of mass motion of a pair from its internal de- 
gree of freedom. The former is described by the standard 
two-particle T— matrix. The imaginary part of the latter 
is essentially the spectral function of a pair with momen- 
tum P, but it does not give any information about the 
internal structure of the pair. The dependence of the pair 
propagator on the internal distance between the partners 
of the pair shows up in what we call the " pair structure 
function" (7p(P,p), which gives essentially the probabil- 
ity for the two partners of a pair, with a given center of 
mass momentum P, of being at the distance p from each 
other. We present and discuss our numerical results in 
section HI, where we show the influence of various pa- 
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ramctcrs, such as temperature, chemical potential and 
interaction strength, on gp .The main conclusions from 
our work are summarized in section IV. In Appendix A 
we discuss the form of the one-electron spectral functions 
in the pseudogap region and in Appendix B we relate our 
pair structure function to the total number of preformed 
pairs which is usually introduced in the framework of 
the Landau- Ginzburg description of the superconducting 
phase transition. 



II. THEORETICAL DESCRIPTION OF THE 
PAIRS 

We consider an extended attractive two-dimensional 
Hubbard model given by the Hamiltonian: 

H=-t ^ {cl^cj^ + h.c.) -^U{i-j) rmrij^ 
-/x^K^-hn^i) (1) 



where t is the hopping parameter between nearest neigh- 
bor sites and U{i — j) > is an attractive interaction 
coupling the electronic density riia — cl^Cia at two dif- 
ferent sites i and j. All our calculations are performed in 
the grand-canonical ensemble, using the chemical poten- 
tial /i, rather than the electronic density, as a parameter. 

The Hubbard model has been analysed by many au- 
thors, using various kinds of approximation techniques, 
in order to describe the instability leading to supercon- 
ductivity and to investigate the influence of the attrac- 
tive interaction on the electronic properties, both above 
and below the superconducting transition temperature 
Tc- By increasing the strength of the attraction, it is 
possible to follow the evolution of a variety of physical 
quantities in a cross-over between a BCS-like behaviour 
to a scenario where superconductivity can be related to 
the condensation of preformed pairs.* Our present objec- 
tive consists in analysing, in the framework of (1), the 
structure of the preformed pairs in the normal state for 
T > Tc, in the domain where simulations show the for- 
mation of a pseudogap. One of the current interpreta- 
tions of this unusual feature in the cuprates precisely in- 
vokes the presence of preformed electronic pairs, existing 
up to a temperature T* , that can be substantially higher 
than Tc- The question of the structure of the pairs below 
Tc will be treated elsewhere. 

Let us consider the operator 



(2) 



defined on the lattice sites i and j with coordinates Xi 

and Xj, which creates a pair of electrons with opposite 
spins. We define the center of mass coordinate of the 
pair and the relative distance , as: 



Pij . 

Although and its adjoint do not satisfy the canoni- 
cal commutation relations for independent particles, they 
can be used in order to describe electronic pair correla- 
tions through the dynamic correlation function: 

x{'rij,Pij,t;rrs,Prs,t') = {-i) (3) 
X {T [Q{rij,pij,t)Q^rrs,Prs,t')] ) , 

where T denotes the time ordering operator. Manifestly, 
X is a two-electron Green's function G2 

G2(l,2;l',2') = -(T[c(l)c(2)ct(2')ct(l')]>, (4) 

with a special choice of arguments (the numbers i = 

1, 2, 1', 2' are a shorthand notation for the generalized 
coordinates Xi, ti, ai of space, time and spin variables): 

Xi^ij ) Pij ) ^rs ) Prs ) O ~ * ^2 (''i j t ; Xj , t ; Xr , t ; Xg , t ) . 

(5) 

A simple and convenient approximation for evaluating 
G2 for a system described by the Hamiltonian (1) consists 
in solving the Bethe-Salpeter equation^: 

G2(1,2;1',2')=G(1,1')G(2,2') - G(l, 2')G(2, l') 

+ i G(l, 1)G(2, 2)[/(l - 2)G2(1, 2; 1, 2) (6) 

taking into account only repeated scattering events be- 
tween two electrons (thus restricting our analysis to low 
density systems) and expressing G2 in term of the "full" 
one-electron Green's function G and of the attraction po- 
tential U. Here and in the following we use imaginary 
time arguments, Q < t < i(3, P = l/iksT), with the 
corresponding time-ordering. We also adopted in (6) the 
convention for which the bar under two identical argu- 
ments denotes integration over space and time. 

Introducing the special choice of variables needed for 
the pair correlation (5) into the general Bethe-Salpeter 
equation (6) and denoting the dynamic correlation func- 
tion (3) by Xiyij - J^rs,Pij,Prs,t - t') EE x{TCi,Pj,Pl,t) 

(where now = ia and Pj = jn are discrete vectors, 
that represent the distance among two sites on a (square) 
lattice of constant a), one obtains 

x{^i,Pj,Pi,t) = -iGirz,Pj - pi,t) + i J dt'J2J2 

m k 

X g{Vi - Trn.Pj - Tfc,* - t')U{Tk) 

X X{^m,rk,Pl,t') (7) 

where Q is the particle-particle bubble: 

g{ruPj,t) = Giri + ^,t)G{ri-^,t), (8) 

describing the motion of two independent particles whose 
properties have been renormalized by the interaction. 
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The spin indices have been omitted, since we are inter- 
ested in the particular configuration of opposite spin of 
the partners of the pairs, and the interaction is indepen- 
dent of the spin variable, so that no spin-flip processes 
have to be considered. 

We now apply a Fourier transformation with respect to 
the variable r, - corresponding to the difference between 
the center of mass coordinates in x - to equation (7), 
introducing the center of mass momentum P, and we go 
from imaginary time to (bosonic) Matsubara frequencies 
Zj,. This yields 

X{^,Pj,Pk,Zu) = -ig{P,Pj -Pk,z^) + i ^ U{ri) 

i 

X G{'P,Pj -ri,z^)x{'P,ri,Pk,Z:,). (9) 

For a local attraction U{r) = UoS{r) equation (9) be- 
comes algebraic and x can immediately be related to the 
particle-particle bubble function Q (see below). When 
the interaction has a finite spatial extension, which is 
needed for treating d-wave pairing, one has either to solve 
the integral equation (9) or to introduce a separable form 
for the Fourier transform of U : 

= ^C/,<^,(k)</.*(k), (10) 

X 

(px being the various characteristic shape functions of the 
potential (the limit of the local potential can be retrieved 
by taking only the mode = 1). In this more general 
case, equation (9) is replaced by 

X(P,p,p',^.) = e-*''(''-''')^(P,k;0,) +i^C/A 

k A 

x5^e-*-VA(k)a(P,k;^.) 

k 

x^e-^'^''''FA(P,k;z.), (11) 

k' 

where Q{P ,'k.; z^,) is the Fourier transform of the particle- 
particle bubble (8) with respect to the the center of mass 
position r, internal distance p and time, and Fx is an 
intermediate-step quantity obtained from the convolution 
of the Fourier transform of the dynamic correlation func- 
tion X with respect to the variables Pi and Pj, and the 
characteristic shape function of the potential: 

F;,(P,k;z.) = ^</.A(k') ^ e^'^^r^ "' 

k' i j 

X x{'P,Pi,Pj,z>^)- (12) 

In the following, we will select one value of A and treat 
two cases separately: 

• s-wave pairing: ^^(k) = 1; 



• d- wave pairing: = cosfc^ — cos ky, 

where we have assumed the lattice constant to be a = 1. 
We will finally aim at the following quantities: 

1. the usual two-particles T-matrix T(P,Zi,), which 
contains all the details of the interaction but does 
not show any dependence on the internal distance 
Pi between the partners of the "pair" . This quan- 
tity is given by 



^ ''"^ 1 -iC/AEkl'^A(k)Pa(P,k;^,)' 

(13) 

2. the spectral function $(P,p, Q) of a pair with cen- 
ter of mass momentum P, energy fl and internal 
distance p: 

$(P,/t),f}) = Im{x{P,p,p,n + iO+)}; (14) 

3. the "pair structure function" gp(P,p), that gives 
the internal structure of a pair having momentum 
P and internal distance p: 

gp{p,p) = J dnNB{n)^p,p,n), (15) 

Nb being the Bose-Einstein distribution function. 

The equations (11) and (12) can now be solved in two 
steps. First, taking the double Fourier transform of (11) 
one can relate F to the Fourier transform of the "bubble" 
g with respect to the variable p: 

Fx{P,\^;z.) = Mk f^^;'"^ giP,k;z.). (16) 

Inserting this result into (11) yields: 

x{P,p,p,z,) = -ig{P,0,z,) +T{P,z,) 

xg{p,p,z^)g{p,-p,z,), (17) 

with a "bubble" weighted by the form function (px- 

g{P,p, z,) = e*-''</)A(k)a(P, k, z,) . (18) 

k 

The two-fold effect of the interaction potential is 
clearly visible in this result. On one hand, it leads to 
a "renormalized" one-electron propagator G and thus to 
a bubble g describing the motion of two independent 
"renormalized" quasi-particles. On the other hand, the 
T-matrix represents the scattering processes that finally 
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lead to the formation of the preformed pairs. Follow- 
ing expression (15) wc can thus construct a function, gp, 
that represents the internal structure of two independent 
particles with center of mass momentum P and which is 
given by the first term in (17): 

gUP,p) = J d07VB(0)/m{x°(P,/£),/£),^.)U.^n+io+} 

= J dnNB{n) Im{g(P,0, z^)U^^a+^o+} ■ (19) 

Being independent of p, this contribution forms a con- 
stant "background" to the structure of gp, resulting from 
the direct two-particles scattering. For this reason, wc 
subtract this background contribution from gp, so that 
instead of the defintion (14) we will use the "rescaled": 

$(P,p,n) = Im{x{P,p,p,n + ir,) - x°(P,o,^ + iv)}- 

(20) 

However, one has to bear in mind that possible nega- 
tive values of gp simply mean a reduction of the number 
of corresponding pairs with respect to this independent 
particle background. 

We do not attempt a fully self-consistent calculation of 
G - and thus of 5 - and %• Rather, we approximate the 
one-electron spectral function in the pseudogap regime 
we are interested in by a BCS-like form with a finite 
linewidth 7. The corresponding result for Q is derived in 
Appendix A. 

In the following section we will discuss the general be- 
havior and present numerical results for the T— matrix 
and the pair structure function for the s— and d— wave 
pairing symmetries. 

III. NUMERIAL RESULTS 

Our numerical calculations are based on the tight- 
binding spectrum resulting from the kinetic energy part 
of the Hamiltonian in (1) with a nearest neighbor hop- 
ping t in two dimensions. We use the total bandwidth 
W = 8t as our energy unit. The zero of the energy is in 
the centre of the band. Assuming that the underdoped 
cuprates with a pronounced pseudogap are situated in 
an intermediate coupling regime of the attractive Hub- 
bard model^" we chose an interaction strength U is of 
the order of 0.5 W. Numerical simulations^^ of the at- 
tractive Hubbard model with a local interaction of this 
strength in two dimensions point to a critical temper- 
ature Tc on the order of 0.01 W for s-wave symmetry. 
This result has to be taken as a qualitative indication, 
since it is not clear whether quantum Monte Carlo cal- 
culations can really represent correctly the Kosterlitz- 
Thouless-type (KT) phase transition that is expected in 
two dimensions. Prom Quantum Monte Carlo calcula- 
tions, the gap value in the s-wave case for the interaction 
strength half of the bandwidth, is of the order of t/2. 



Based on these numerical results from simulations, 
we choose a gap parameter value for which the critical 
temperature is approximately 0.2A in the s— wave case. 
In the d— wave case, however, the interaction strength 
should be of order of the bandwidth to obtain reasonable 
values for the critical temperature,^^ and still maintain 
the proportion between the critical temperature and the 
value of the pseudogap approximately equal to 0.2. 

In order to incorporate the effect of the finite lifetime 
of the pair excitations without washing out the effect of 
the pseudogap, we choose the finite linewidth 7 (see Ap- 
pendix A) to be of the order of A/2. In this way, in fact, 
one can still see rather sharp features in the pseudogap 
region, while for bigger values of 7 they are strongly re- 
duced. 

The relevant temperature interval of the pseudogap 
regime is determined in the framework of our description 
of pairing. In the usual T-matrix approximation, which 
results from the Bethe-Salpeter equation (6), the super- 
conducting instability occurs at the temperature where 
the "Thoulcss criterion" is fulfilled: 

r-i(P = 0, = 0) = 0, (21) 

i.e., where pairs with zero centre of mass momentum be- 
gin to form a true bound state. This is a "mean-field-like" 
criterion that corresponds to our approach to pair struc- 
ture, which is, of course, less precise than KT, but the 
latter would require to go more deeply into the analysis 
of phase fiuctuations. 

In Figure 1 we show the effect of temperature on the 
real and the imaginary parts of the T-matrix r(0, fl + 
iO~^) for the case P = 0. As the temperature approaches 
the transition value, the peaks in the real and imaginary 
part of the T— matrix become sharper, and move toward 
LU ~ 0. The inverse of the width of the imaginary part 
peak can be interpreted as the lifetime of the preformed 
pairs. We obtain for the lifetime the values 25, 5.5, and 
1.8, for T/Tc equals to 1.5, 3 and 6, respectively. The 
lifetime is expressed in unit of the inverse of the Fermi 
energy Ep, where Ep = fj, - Ebottom and Ebottom is the 
bottom of the band energy value. 

For our calculation, we will then consider the two tem- 
perature ratio: T/Tc = 1.5 and T/Tc = 3 to constrast 
the cases in which the quasiparticle features are more or 
less pronounced. 

We study the behaviour of the function gp for vari- 
uos values of the pair internal momentum. Moreover, we 
also investigate the influence of temperature, chemical 
potential and interaction strength on the shape oi gp. 

For comparison, we will also show the results for A = 
for s— and d— wave interaction symmetries. Although the 
pair structure function is defined on a lattice, we present 
it as a continuous function, as a guideline for the eye to 
follow its dependence on the relative distance. 
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A. s— wave pairing 

We begin by considering the simplest case of A = 0, in 
order to show the general tcndcncc of the pair-structure 
function without the complication of the influence of the 
pseudogap on the pairs. In Figure 2, we show a typi- 
cal set of curves for s— wave symmetry pairing for the 
case A = 0. The function is calculated for the relative 
distance directed along the x axis in the Brillouin zone. 
In the s— wave case, however, the pair structure func- 
tion is isotropic, reflecting the spatial symmetry of the 
order parameter, so that gp(P,px) is representative of 
the behaviour along the other directions. Moreover, the 
pair structure function depends mainly on the modulus 
of the centre of mass momentum, while the variations of 
the direction of P have no relevance for gp. As can be 
easily seen, the effect of increasing pair internal momen- 
tum is to reduce the number of pairs and to shrink their 
size (which can be seen as a sort of "relativistic length 
contraction"). It is apparent that the extension of the 
pairs in all the cases here considered is of the order of 
the lattice constant. The effect of the temperature on 
the distribution function is seen by comparing Figure 2a 
and 2b. As expected, raising the temperature reduces 
the number of pairs. A slight decrease of the avarage ex- 
tension of the pair can also be observed. We notice also 
that increasing the temperature tends to wash out the 
difference of population of pairs with different internal 
momentum. In fact, due to the increase of the kinetic 
energy for particle, it becames easier to form pairs with 
finite center of mass momentum, so that the number of 
the pairs with P 7^ increases, even if the net effect of 
increasing temperature is to lower the total number of 
pairs. 

In Figure 2c we show gp for a different value of the 
chemical potential /x. The comparison is made with Fig- 
ure 2a, at a fix temperature ratio with respect to the crit- 
ical one: T/Tc = 1.5. As the chemical potential moves up 
from the bottom of the band (i.e. for increasing density, 
but still in the limit of validity of our T— matrix approx- 
imation), the number of pairs increases, while the size of 
the pair slightly decreases. Fixing the chemical potential 
and for a stronger interaction, we can observe that the 
number of pairs increases, while their spatial extension 
does not vary appreciably (Figure 2d). 

We present now the results for the pseudogap region, 
i.e. for A ^ (Figure 3). The effect of a pseudogap 
opening in the single particle spectrum is to increase the 
number of pairs and to enlarge their size for zero internal 
momentum, while for finite value of P it seems that the 
opening of a pseudogap suppress the pairing, as one can 
see comparing Figures 2a and 3a. As already mentioned, 
we choose the parameters for the pseudogap region from 
the results of Monte Carlo simulation. ^° The dependence 
of the (pseudo)gap parameter A on the chemical poten- 
tial and on the interaction has been stated by fixing the 
relative magnitude of the critical temperature and the 



pseudogap amplitude to be the one tipical for cuprates. 

As wc noticed, the opening of a pseudogap strongly de- 
presses the number of pairs having a finite momentum. 
This behaviour becomes even more evident when the 
strength of the interaction increases (Figure 3d), which in 
turn implies a bigger value of the pseudogap. A part from 
the different behaviour at finite momentum, the qualita- 
tive dependence of the pair distribution function on the 
parameters in presence of a pseudogap does not differ 
from the case A = 0. Also, we have not found any sub- 
stantial variation of the avarage internal distance of the 
pairs by varying the interaction strength. This is in part 
due to the fact that we confine our study to the pseudo- 
gap region, so that the range of variation for the inter- 
action strength is not sufficient to crossover from Bose- 
Einstein condensation to BCS behaviour. In Figure (4) 
we show the variation of the pair internal distance with 
respect to variation of the interaction strength for the 
case A = 0, where we have normalized the pair struc- 
ture function to its zero internal distance value, as we 
want to emphasize the trend of the functions. As it can 
be easily seen, the internal distance has doubled as U has 
passed from 0.61^ to 0.2W. 

This evolution is similar to what has been obtained by 
Paredes and Cirac^^ for a system of bosonic atoms in an 
optical lattice, which is mapped into a fermionic system 
via a Jordan- Wigner transformation with a local attrac- 
tion V. The pairs are formed according to a variational 
wave function, whose spatial extension changes by chang- 
ing the attraction strength, and results in localized pairs 
for large value of V , and delocalized pairs for small V. 

In the context of cuprate supercondutors, it is inter- 
esting also to consider the extended s— wave symme- 
try, which has been recently proposed to explain some 
experiments.^^ In Figure 5 we present the pair struc- 
ture function calculated in the pseudogap region for 
0A(k) = coskx + cos ky. We can see that in this case 
the maximum of the distribution is shifted from zero in- 
ternal distance to a finite distance of the order of one 
lattice constant. However, the effects of the variation 
of the parameters on the pair structure function are the 
same as in the full gapped s— wave. 

B. d— wave pairing 

The d— wave case does not differ too much in the qual- 
itative behaviour from the previous cases. Nevertheless, 
in the d— wave case, gp vanishes for zero internal dis- 
tance, as it should be from the symmetry of the interac- 
tion and of the pseudogap parameter, and has its max- 
imum at a distance of the order of a lattice constant, 
slightly higher than the extended s— wave case. The fact 
that the pairs have again typically the size of a lattice 
constant is a consequence of our nearest neighbor type 
of attraction. As for the s— wave case, we first analize 
the case A = (Figure 6). Besides the general features 
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alrcaady mentioned, we can see that a secondary peak 
appears in the zero momentum distribution function at 
a distance of three or four lattice constant. Moving pairs 
show a shght contraction and they are less numerous than 
pairs at rest. As the temperature is raised, the distribu- 
tion function is suppressed for P = and enhanced for 
P ^ 0, as in the s— wave case (compare Figure 6a and 
6b). 

By lowering the interaction strength (Figure 6a and 6d), 
the number of pairs with internal momentum P tends to 
diminish strongly. Moving the chemical potential more 
inside the band (Figure 6a and 6c) has as a result an 
enhancement of the number of pairs. In this case, the 
increase in the population of pairs with finite momentum 
is remarkable, in contrast to the behaviour obtained when 
the interaction strength is changed, when no noticeable 
change in the relative populations for different P appears. 
The result obtained by varying the chemical potential can 
be attributed to a greater avarage kinetic energy due to 
the spatial confinement of the particle at higher densities. 
The form of gp shows the typical Friedel oscillations, the 
wave length of which depends on electronic density (as 
can be seen when changing the value of /i). 

The opening of a pscudogap (Figure 7a) strongly en- 
hances the number of pairs for total momentum P — 0, 
while suppresses the number of pairs with P ^ 0. No 
change in the pair size is instead noticeable. Another 
fact to be noticed is that when a gap opens we need to 
move toward stronger interaction value in order to have 
a reasonable value for the transition temperature. This 
is due to the fact that a finite value of the interaction in 
the 2D Hubbard model is needed for the d— wave case to 
have a bound state in the related two-body problem. 
Moreover, the Friedel oscillations are in this case hidden 
by the oscillation of the d— wave symmetry of pseiidogap 
parameter, of typical period n/a. 

As expected, we obtained that the d-wave pairs are 
elongated in the direction of the coordinate axis, whereas 
their width shrinks to zero for p along the diagonal 
{Px = Py) (see Figure 8). As in the s— wave case, the 
pair structure function is not very sensitive to the direc- 
tion of P, but depends only on its magnitude. 

The effects of temperature and interaction strength on 
the pair structure function are not different in the general 
behaviour from the case with no pseudogap. To increase 
the chemical potential depresses strongly the zero mo- 
mentum pair distribution. On the contrary, the popula- 
tions of pairs with finite momentum increase. To explain 
this enhancement, we can invoke the same explanation as 
in the A = case (namely, the increase of kinetic energy 
due to confinement). We can also note a decrease of the 
avarage extension of the pair, as the secondary peak in 
the pair structure function in figure 7c has almost com- 
pletely disappeared. It is to note, however, that the total 
number of pairs seems to decrease with increasing chemi- 
cal momentum, despite the fact that the total number of 
particles is higher. This could be attributed to a break- 
down of our T— matrix approximation, or to the fact that 



we are approaching a region in which the pseudogap is 
closing. 

IV. DISCUSSION 

In this work, we have investigated the internal struc- 
ture of preformed pairs in the pseudogap region of high 
temperature superconductors. We assume that this re- 
gion can be described by an attractive Hubbard model. 
In this framework, we have studied the Green function 
for the operator which creates a pair of particles whose 
partners are at a certain relative distance. We applied 
the Bethe-Salpeter approach and introduce a T— matrix 
formalism, which takes into account the relevant inter- 
action channel of repeated two-particle scattering. Our 
calculations show that in a system described by an attrac- 
tive Hubbard model, the opening of a pseudogap in the 
one particle spectrum increases the probability of find- 
ing two particles at a distance of the order of the lat- 
tice spacing, with respect to the case in which no gap 
is present. We have also found that in the range of val- 
ues that we have considered for the parameters (and that 
should correspond to the pseudogap region) the average 
extension of the pair is of the order of a lattice constant. 
The interparticle distance docs not change appreciably 
by varying the interaction strength, but this is not too 
surprising since we have not investigated the BCS limit 
of our model, where the pseudogap region is believed not 
to be present. Nevertheless, in the broken symmetry 
phase the dependence of £,pair at T = on the inter- 
action strength is instead quite strong^'^'', such that one 
could presume that the anomalous Green functions might 
in fact bear the major effect to the extension of the pair. 
The effect of increasing temperature is to reduce the to- 
tal number of preformed pairs and to shorten their in- 
ternal distance. Also, by increasing the avarage kinetic 
energy of the particle, it changes the relative weight of 
the populations corresponding to different centre of mass 
momentum P, depressing the one with P = more than 
the ones with P 7^ 0. The influence of the chemical po- 
tential on the pair distribution function reflects the influ- 
ence of the density of particles. It is not surprising, then, 
that moving the chemical potential more inside the band 
(i.e., for higher densities) increases the number of pairs 
with higher momentum: the kinetic energy of the particle 
should in fact increase as well, and so does the avarage 
momentum of the pairs, since the number of pairs for a 
given momentum is given by the integral over the inter- 
nal distance of the pair structure function (see Appendix 
B). 

One interesting question concerns the connection be- 
tween the preformed pairs and the structure of Cooper 
pairs below Tg. We in fact expect to recover in the broken 
symmetry phase the usual BCS to Bose-Einstein conden- 
sation behaviour, which obviously affects also the average 
internal pair distance. 
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Another interesting observation related to the pair cor- 
relation function, which is equivalent to our pair struc- 
ture function integrated over all the centre of mass mo- 
menta, is the fact that in the broken symmetry phase it 
is connected to the correlation energy.^* This energy, in 
turn, plays a key role in the thermodynamics of the phase 
transition itself. The real space picture of this pair cor- 
relation function, reported in^^ looks similar to what we 
obtained for gp{P,p) above Tc both for s— and d— wave 
case as can be seen by comparing with Figure 8 (we point 
out that in contrast to reference^^, we reported the val- 
ues of the pair structure function not only on the lattice 
sites, but also on intermediate positions). 

Besides this aspect related to superconductors, another 
possible relevance for the internal structure of Cooper 
pairs is in the very active field of systems of ultracold 
atoms. In fact, while it is very difficult to measure the in- 
ternal structure of the Cooper pairs in a superconductor, 
due to the very complicated structure of the cuprates, in 
the case of atoms trapped in optical lattice it is possible 
to carefully control all the parameters for the interac- 
tion among particles^^. Moreover, it should be possible 
to actually measure the extension of the Cooper pairs, 
for example, by analyzing the response to a Raman laser 
which couples to the internal states of the pair, as sug- 
gested in^^. 

In order to make a more meaningful comparison with 
the cited papers dealing with the superconducting phase, 
we should however calculate the pair structure function 
below Tc. That will be the subject of a successive work. 

We thank S. Sharapov, P. Fieri and G. C. Strinati 
for helpful discussions. This work was supported by the 
National Competence Centre "MaNEP" of the Swiss Na- 
tional Science Foundation. 



APPENDIX A: SPECTRAL FUNCTION IN THE 
PSEUDOGAP REGION 

To mimic the presence of a pseudogap, we consider the 
form of the BCS Green function with a finite linewidth. 
Due to the absence of long range order the "anomalous" 
BCS Green function is zero. We have then: 



G(k, 



1 



z - ^(k) + iV ■ 



A^(k) 
z-(,(\^)+iT 



(Al) 



where ^(k) = —2t{coskx + cosky) is the single parti- 
cle disperion, l/F is the excitation lifetime and A(k) is 
the pseudogap. It is possible to rewrite (Al) in terms 
of the quasiparticles excitation energy, that is i?(k) = 



z + e(k) + iT 



{z + E(k)+iT){z-E{k)+ir) 



z-E{k) + ir z + E{k) + ir' 



(A2) 



where = (^(k) + ^(k))/(2i;(k)) and vl = {E{k) - 
^(k))/(2iJ(k)). From (A2) we can deduce the form of 
the spectral function: 



A{k,z) = - 



+ 



TT V(-2--E(k))2-Kr2 {z + E{-k))'^ 

(A3) 

The particle-particle bubble can be written in term of 
the spectral function of the single particle Green func- 
tions: 

g{P, Kz)=J2 G(k + |, ^.) G(k - |, z - z,) 
Z — UJ — to' 

Replacing the expression for the single particle spectral 
function (A3) into the previous equation and performing 
the substitution z — > Q.-\-i5, we can obtain the expression 
for the spectral function of the particle-particle bubble 
/(P,k,f^): 



/(P, k, n) = ImGiP, + iS) = 
2 / uju^r 



n \{n - {El + E2))^ + {n+{Ei+E2))^ +4:T2 

+ (fi _ {El - E2)f + 4r2 + (fi + (^1 - E2)f + 4r2 

(A5) 

where we used the shorthand notation 1 = k -|- P/2 and 
2 = k — P/2. The parameter that enters our calculations 

and that mimics the pseudogap is then the linewidth 7 = 
2r for the excitation described by the particle-particle 
bubble. 



APPENDIX B: NUMBER OF PREFORMED 
PAIRS 

In this Appendix we relate the definition of the pair 
correlation function (3) to the number A^c.p. of fluctuating 
Cooper pairs above Tc- If the analysis is carried out 
close enough to the critical temperature, one can define 
the density of fluctuating Cooper pairs starting from the 
Ginzburg-Landau (GL) theory^^. Above Tc, the order 
parameter ip{p) has a fluctuating origin and depends on 
spatial coordinates even in the absence of magnetic field. 
The fluctuating thermodynamical potential is then: 



dr 



e I V'(r) P +^ I Hr) 



2a 



+ Vd I V'^i'(r) 



(Bl) 



where e = In (T/Tc) <C 1, a = l/{4mrj£)) and 77/3 is 
the phenomenological GL constant. The GL order pa- 
rameter undergoes equilibrium fluctuations which, close 
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enough to the transition temperature, can be described 
by a probabihty distribution^^: 



V oc exp = exp 



(B2) 

where Fourier transform of the order parameter with re- 
spect to the spatial variables has been carried out. The 
average value for the square of the fluctuating GL order 
parameter is then: 



/nk)d|^kP 



(B3) 



The number of fluctuating Cooper pairs Nc.p. is deter- 
mined by the average value of the square of the GL order 
parameter modulus^": 



Nc.p. = 5^(|Vri^>expi(k.r) 



r->0 



(B4) 



The number of fluctuating Cooper pairs has been ob- 
tained measuring the Knight shift and the nuclear 
spin relaxation under magnetic fleld by Zheng and co- 
workers^^. These quantities reflect indeed the pseudogap 
behaviour, and their reduction above Tc is strongly mag- 
netic fleld dependent accordingly to a scaling relation 
that depends on the effects of the Cooper pair density 
fluctuations^"'^. 

To relate the deflnition (B4) to our pair correlation 
function (3), we consider the partition function for the 
Hamiltonian (1). For the sake of simplicity, we consider 
the case of on-site interaction U{i — j) = U 6i j, where S 
is the Kronecker delta function. By means of a Hubbard- 
Stratonovic transformation we can rewrite the partition 
function in term of the classical fleld A(P, r) and its com- 
plex conjugate: 



Tre 



-fSH 



Tr e-/^"" T j VA VA* exp 

J2\U\-'\A{P,r)f + (gt(p,^)A(P,r)+ 
p 

go(P,r)A*(P,r))]dT, 



Jo 



(B5) 



where T is the time-ordering operator and Qo(P,t) is 
the Fourier transform with respect to spatial coordinates 
and time of the local pairing operator defined in (2) for 
i = j. Hq is the non interacting part of the Hamiltonian 
H (Eq. 1). 

In order to obtain the two-particle Green function we 
are interested in, we introduce the source fleld, ry(P,r) 
and its complex conjugate, which couples to the pairing 
operator, to obtain the action term (replacing the expo- 
nent of (B5)): 

S= [ ^''^|f/|-MA(P,r)p + Qj(P,r)(A(P,r)+ 
Jo p 

+ 77(P, r)) + Qo(P, r) (A*(P, r) + r)) dr . (B6) 



If we now differenciate the partition function Z with re- 
spect to the flelds t] and rj* , we obtain the Green function 
for the local pairing operator: 

X 



6'n*(P',T') (5r?(P,r) 



(T[gS(P,r) Qo(P',t')]). (B7) 



We can now define a new field: ^(P,t) = A(P,t) -|- 
?7(P,r), so that Eq. (B6) becomes: 



Si 



Jo D 



\Ur\aP,T)-v{P,T)\' 



+ QliP, r) C(P, r) + Qo(P, t)^ (P, t) dr . (B8) 

Differentiating Z again with respect to rj and r]* we ob- 
tain the equality: 



(T \Ql{P,T)QoiP',T')\) \r'=r 



(A*(P,t)A(P',t)) 



C/2 



I 

W\ 



(B9) 



where we have set r' = r in order to discuss the relation 
among this expression and the equation for the number 
of Cooper pairs derived from the (time independent) GL 
theory. Above the critical temperature, the fluctuating 
GL order parameter is proportional to the classical field 
A(P, r) we have introduced. Equation (B9) essentially 
confirms formula (B4), which represents the number of 
fluctuating Cooper pairs that have a fixed distance be- 
tween the components, namely p = for an on-site in- 
teraction. 

We want now calculate the Green function for the Q 
operator for an internal distance p. To this end, we must 
insert in (B5) a source fleld that depends on p, namely 
?7(P, P,t): 

^0 p 

+ Qo(P,r)A*(P,T) V (Qt(p,p,^)^(p,p,^) 



+ Q(P,p,t)v*(P,P,t)) dr. 



(BIO) 



The correlation function x(P,p,p',t) can be obtained by 
differentiating with respect to the source flelds: 

X X 

TogZ|^=0,r;*=0 = 



5r]*{P',p',T')S7jiP,p,T) 

= {T[QHP,p,t) Q(P',p',t')]). (Bll) 

In order to obtain the relation between (Bll) and the 
classical fleld A, we perform a second order cumulant 

expansion for the trace over electronic degrees of fredoom 
in Z. The odd power terms of the expansion vanish. We 
then obtain: 



8 



-I 



Z = Zo DADA* exp 



^\UnA{P,r)fdT 



X exp 
where: 



C{P,T;P',T')dTdT' 



p,p' 



(B12) 



C(P,T;P',r') = 

^(T [Q\P,p,T)Q{P',p',T')])MP,P,rW{P',p',r') + 



p,p' 



+ ^(T [Q\P,p,T)Qo{P',T')]},r,{P,p,T)A*{P',T') + 
p 

Y,{T [0S(P,T)Q(P',p',r')])cA(P,T)7?*(P',p',T') + 



+ 



p' 



+ {T [Ql{P,T)QoiP',r')\UA{P,T)A*{P',T') + h.c. . 

(B13) 

{■ ■ ■)c is the second order cunmlant, which in our case 
coincides with the avarage over Hq : (•••)o- Applying 
the definition (Bll) to the expression (B12), we obtain: 

{T[Q''{P,p,t)Q{P',p',t')])2 = 

-{T[QHP,p,r)Q{P',p',T')])„+ f fdT"dr"' ^ 



X (T [gt(P", 0, r")g(P',p', r')] )o(T [Qt(P, p, r) 

X g(P"',0,T"')])o(A(P",T")A*(P"',T"'))2. (B14) 

The trace (• • • )2 is taken over the second order cumulant 
expansion of the total Hamiltonian. 

We specialize the equation (B14) to the case P = P' 
and p — p' , and apply Wick's theorem to the trace over 
Hq. Recalling the definition of the dynamic correlation 
function % (Eq. 3) and of the particle-particle bubble Q 
(Eq. 8), we can make the following identification: 

(T [Qt (p, r)Q(P, p, r')] )2 = i x(P, P,t - r') (B15) 

(T [Q\P, p, t)Q{P,p, t')] )o = - g(P, 0, r - r') (B16) 

(T[Qt(p,o,r)Q(P,p,r')])o = - g{P,p,T - r') . (B17) 

Performing a Fourier transformation with respect to 

time, we finally recover from equation (B14) the result 
obtained with a Bethe-Salpeter approach; 

= -^j dre'-^-(r[Qt(P,p,r)Q(P,p,0)])2 

= -ia(p,o,^,) + g{p,p,z,,)g{p,-p,z,) 



u 



1 - iC^Ek e(P,k;^.) 



(B18) 



where the Fourier transform of the trace of the classical 

fields calculated with respect to the second order cumu- 
lant expansion distribution has been identified with the 
T— matrix: 



j dTe''^^A{P,T)A*{P,0))2=iT{P,z^) 



iU 



1 - i?7Ek e?(P,k;z,) 
(B19) 



Equation (BIB) is equivalent to (17) for the case of a local 
potential (i.e. (pxCk) = 1). 

Expression (B18) calculated for p = is equivalent to 
the r.h.s. of equation (B9). This shows again that the 
number of Cooper pairs calculated via eq. (B4) gives the 
number of Cooper pairs with fixed internal distance, i.e. 
zero internal distance for an on-site potential. 
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FIG. 1. Imaginary part (left panel) and real part (right 
panel) of the T-matrix for the d— wave case for different tem- 
peratures t = T/Tc and P = 0. 





FIG. 2. "Pair structure function" gp{p) for p = (p, 0) (in 
units of a) and A = and: (a) T = 1.5Te, fi = -0.4H^, 
U = 0.5W; (b) T = 3Tc, n = -OAW, U = 0.51^; (c) 
T = 1.5Tc, fj, = -0.3W, U = 0.5W; (d) T = 1.5T^, 
fx = -OAW, U = 0.6W. 



10 



(a) 

Pj,'-0 
P^=7i:/ 4 
P^=7i:/2 

P„=7C 




?^=7C/^4 

P^=7C/2 
?„-37t/4 



0.5 1 1.5 2 2.5 3 




FIG. 3. "Pair structure function" gp{p) for: (a) 



A = 0.04W^ T = 1.5Tc, ^ = -OAW, U = 0.5W": 
A = 0.04^-^^, T = STc, /i = -0.41^, [/ = 0.5W 
A = 0.062VK, T = 1.5r„ ^ = -0.3H/, U = 0.5W 
A = 0.067^^, T = l.5Tc, n = -OAW, U = 0.6W. 
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FIG. 5. "Pair structure function" gp{p) for: 

A = 0M7W, T = l.5Tc, fi = -OAW, U = 0.5W 

A = 0M7W, T = 3Tc, /j, = -OAW, U = 0.5W 

A = 0.054W^, T = 1.5Tc, = -0.3W^, U = 0.5W 
A = 0.078Ty, T = 1.5Tc, m = -0.41^, U = 0.7W. 
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FIG. 6. "Pair structure function" gp{p) for d— wave case 
and A = 0: (a) T = 1.5Te, U = O.&bW, /x = -OAW; 
(b) T = 2,Tc, U = 0.65W", fi = -OAW; (c) T = 1.5Tc, 
f/ = 0.65Vy, n = -0.3W; (d) T = 1.5T,, (7 = 0.63Vy, 
/u = -0.4W. 
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FIG. 7. "Pair structure function" gp{p) for: (a) 

A = 0.05W, T = l.STc, M = -0.4Ty, [/ = 0.7W; (b) 
A = 0.05W^, T = STc, /i = -OAW, U = 0.7Vy; (c) A = Q.IW, 
T = 1.5Tc, n = -0.3W, U = Q.7W; (d) A = QMW, 
T = 1.5Tc, n = -OAW, U = 0.68W. 
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FIG. 8. Spatial structure of "pair structure function" 
pp(0,p) for the d-wave symmetry (A = Q.05W, T = 1.5Tc, 
H = -OAW, U = 0.7W) 
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